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^ , ABSTRACT 

Ti^lj- ■ Wc present an analysis of the small-to-intermediate scale clustering of samples of 

Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) and the 
2dF-SDSS LRG and QSO (2SLAQ) survey carefully matched to have the same rest- 
' frame colours and luminosity. We study the spatial two-point auto-correlation function 

\ in both rcdshift-spacc (^(s)) and real-space {£,{r)) of a combined sample of over 10,000 

LRGs, which represent the most massive galaxies in the universe with stellar masses 
> lO^^/i~^M0 and space densities ~ 10~^/i'^Mpc~'^. We find no significant evolution 
in the amplitude (ro) of the correlation function with redshift, but do see a slight 
^ ' decrease in the slope (7) with increasing redshift over 0.19 < z < 0.55 and scales of 

d I 0.32 < r < 32/i~^Mpc. We compare our measurements with the predicted evolution of 

dark matter clustering and use the halo model to interpret our results. We find that 
our clustering measurements are inconsistent (> 99.9% significance) with a passive 
model whereby the LRGs do not merge with one another; a model with a merger rate 
of 7.5±2.3% from z — 0.55 to 2; = 0.19 (i.e. an average rate of 2.4% Gyr^-'^) provides 
a better fit to our observations. Our clustering and number density measurements are 
consistent with the hypothesis that the merged LRGs were originally central galaxies 
in different haloes which, following the merger of these haloes, merged to create a 
single Brightest Cluster Galaxy (BCG). In addition, we show that the small scale 
clustering signal constrains the scatter in halo merger histories. When combined with 
measurements of the luminosity function, our results suggest that this scatter is sub- 
Poisson. While this is a generic prediction of hierarchical models, it has not been tested 
before. 

Key words: surveys - cosmology: observations - cosmologydarge-scale structure of 
Universe - galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: haloes 
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1 INTRODUCTION 

In recent years, the evolution of massive galaxies in the uni- 
verse has received much attention because of the possible 
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tension between observations of the abundance and cluster- 
ing of such galaxies, as a function of redshift, and predic- 
tions from popular hierarchical models of galaxy evolution. 
Naively, in a Cold Dark Matter (CDM) dominated universe, 
one would expect the most massive galaxies to form last 
through the hierarchical merging of smaller galaxies. This 
behaviour is illustrated in the r ecent high resolution sim- 
ulations of iDe Lucia et ah I (I2OO6'), which include the latest 
semi-analytical formalism and account for feedback from ac- 
tive galactic nuclei (AGNs). In these simulations, the stars 
in the most massive galaxies are formed at high redshifts, 
but their stellar mass is only assembled into a single sys- 
tem at relatively late times through "dry mergers", i.e., 
major mergers of gas-poor galaxies with little or no asso- 
ciated star-formation . For example, in Figures 4 and 5 of 
iDe Lucia et~al] l|2006l ). the simulations shows that for low 
redshift elliptical galaxies, with masses > 10^ ^Mq, 80% of 
their stars are formed at a median redshift of 2 ~ 2.5, but 
80% of the stellar mass is only put in place by 2 ~ 0.3. 
Likewise, the simulations show that galaxies with masses 
> 10^^ Mq have multiple large progenitors and cannot be 
formed through a single major merger of two large galaxies. 

These recent AGN-feedba c k mo d els of galax y evo- 
lution (see also [Croton et al.l I2OO6I : iBower et all I2OO6I : 
iHopkins et al.|[2006l ) solve the apparent inconsistency of the 
old ages of stars in massive galaxies (both in and outside 
galaxy clusters) and the late assembly of such galaxies in a 
A-dominated CDM universe. However, they appear to be in 
conflict with recent observations of the luminosity function 
and clusterin g of massive ellipti cals as a function of redshift. 
For example. fWake et al.l l|2006f ) showed that the lack of evo- 
lution of the luminosit y function of Luminou s Red Galaxies 
(LRG s; as defined in lEisenstein et al.l I2001I : ICannon et al.l 
I2OO6I ) put an upper limit on the amount of allowed evolu- 
tion in these massive galaxies, i.e., at least half of the LRGs 
at low redshift (2 ~ 0.2) must already have been well as- 
sembled (with more than half their stellar mass in place) by 
z r-^ 0.6. This is in excellent agree ment with other lum inos- 
ity function studies. For example, iBrown et al.l |2003) used 
data from the NOAO Deep Wide-Field Survey (NDWFS) 
and the Spitzer IRAC Shallow Survey to show that "~ 80% 
of the stellar mass contained within today's 4L* red galaxies 
was already in place at 2 = 0.7". These observational con- 
straints are barely consistent with the semi-analytical CDM 
simulations discussed above. 

The clustering of ma ssive ellipticals provides an addi- 



tional test of the models. iMasiedi et al.l (|2006l ) argue that 
the small-scale clust ering of LRGs fr om the Sloan Digital 
Sky Survey (SDSS; lYork et al.l I2OO0D suggest that LRG- 
LRG mergers (i.e. a major merger of two equally massive 
systems) were not important for the mass growth of LRGs 
below 2 = 0.36. More recentlv. lMasiedi et al] l|2007h used the 
LRG-galaxy cross-correlation function to study the small- 
scale clustering of LRGs at 2 ~ 0.25 and concluded that 
LRGs grow in stellar mass at most by ~ 10% between 
0.1 < z < 1 (or approx imately half the age of the Uni- 
verse) . |WhiteI^^^l] (|2007l ) interpretted the evolution in the 
clustering of luminous red galaxies in the NDWFS using the 
halo model - they argue that a third of all satellite galaxies 
(in a halo) disappear over the redshift range 0.5 < 2 < 0.9. 
Since the satellite fraction in their models is of order 20%, 
only about 7% of the galaxies have merged. However, if these 



mergers increase the stellar mass of th e central object, t hen 
this increase can be 25% or even larger. iBell et ahl (|2006l ) re- 
port rapid evolution in the stellar mass of red galaxies since 
2 ~ 1. This apparent discrepancy is probably due to the 
differences in the luminosity distributions of the samples, as 
it is known in clusters that most of the evolution on the 
so-c alled "red sequence" i s at magnitudes fainter than L* 
(see iDe Lucia et all I2OO6I : IStott et all |2007| . and references 
therein). 

In this paper, we expand our earlier study of the evo- 
lution of the LRG luminosity function (|Wake et al.l |2006| . 
Paper I) to include an investigation of the two-point auto- 
correlation function of these galaxies. The key difference of 
this work to that in the literature is the combination of two 
large samples of LRGs from the SDSS and the 2dF-S DSS 
LRG and QSO (2SLAQ) survey (ICannon et al.ll2006D . As 
in Paper I, we are careful to ensure the colour selection of 
LRGs is consistent between these two surveys, thus allow- 
ing a study of this unique population of massive ellipticals 
across the redshift range of 0.15 < 2 < 0.6. In addition, 
this paper uses the halo model to understand the evolution 
of the clustering of galaxies and constrain the merger rates 
of LR Gs. Although the logic is similar to the IWhite et al.l 
l|2007l ) analysis of NDWFS, our halo model is entirely ana- 
lytic, rather than entirely simulation-based . Our analysis is 
complementary to that of lRoss et al.l (|2007l l who study the 
redshift space correlation function of the 2SLAQ sample, 
binned in pair separation parallel and perpendicular to the 
line of sight, and f it both bias i ng an d cosmological param- 
eters to this data. IRoss et "all (|2007l ) conclude that "LRGs 
have a constant space density and their clustering evolves 
purely under gravity", which is consistent with the results 
of Paper I. Here we wish to test if this conclusion remains 
true under a more precise comparison of the evolution of the 
correlation function of LRGs where we accurately account 
for the changing definition of an LRG with redshift.. 

In Section [21 we describe the SDSS and 2SLAQ data 
used in this paper, while in Section [3] we provide details 
of the sample selection used to ensure a consistent defini- 
tion of an LRG across the two samples. In Section [J] we 
present our measurements of the two-point correlation func- 
tion in both real and redshift -space. Section \S\ presents a 
halo model analysis of our measurements and discuss con- 
straining the merger rate in the halo model framework in 
Section [B] We discuss our findings in the context of recent 
work in Section [7] and conclude in Section [S] Throughout 
this paper, we assume a flat A-dominated cosmology with 

= 0.27, Ho = 70km s"^Mpc"\ and ag = 0.8 unless 
otherwise stated. 



2 DATA 

We present in this paper an analysis of galaxies taken from 
both the Sloan Digital Sky Survey (SDSS) and the 2dF- 
SDSS LRG and QSO (2SLAQ) survey. The SDSS survey 
contains 2 main spectroscopic galaxy data sets: the MAIN 
sample and the LRG sample. The MAIN sample consists of 
all galaxies with a Galactic extinction corrected petrosian r 
magni tude r^et < 17.77; t his results in a median redshift of 
~ 0.1 (jStrauss et al.ll2002D . The LRG sample uses a series of 
colour and magnitude cuts with the aim of selecting LRGs 
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out to 2 ~ 0.5 (see Eisenstein at al. 2001 for details of this 
sample). Here we only consider the Cut I LRG sample, which 
has a magnitude limit rpet < 19.2 and is designed to select a 
pseudo volume-limited sample of LRGs, with Mr ^ —21.8 
and 0.15 < z < 0.35. At low redshift there is considerable 
overlap between the MAIN and LRG samples. We select 
th ese two samples of galaxies from the SDSS Data Release 
5 lAdelman-McCarthy et al.ll2007l '). 

The 2SLAQ LRG survey was designed to extend the 
SDSS LRG sample to z ~ 0.7. The LRGs were again se- 
lected with colour and magnitude cuts using the SDSS imag- 
ing. Spectra were obtained with the 2dF spectrograph on 
the Anglo- Australian Telesco pe. Full details of the selection 
and observations are given in I Cannon" et all (|2006h . The fi- 
nal LRG sample contains over 11000 LRG redshifts, cover- 
ing 180 deg'^ of SDSS imaging data with over 90% of these 
galaxies within the redshift range 0.45 < z < 0.7. The tar- 
geted LRGs were s plit into three subsamples as detailed in 
ICannon et all ll2006il . with the primary sample (Sample 8) 
accounting for two thirds of these. We only focus on Sample 
8 in this paper due to its high completeness and uniform se- 
lection. The overall success rate of obtaining redshifts from 
the 2dF spectra for Sample 8 LRGs is 95%, while the centers 
of the 2dF fields were spaced by 1.2°, resulting in an over- 
all redshift completeness of sample 8 LRG targets of ~75% 

across the whole survey area (Paper I). 

Although the SDSS magnitude system (iFukugita et al 



199€ ) was designed to be on the AB scale (|Oke fc Gunn 



19831 ). the final calibration has differences from the proposed 
values by a few percent. We have applied the corrections 
mAB = msDSS + [-0.036, 0.012, 0.010, 0.028, 0.040] for u, g, 
r, i, z respectively (Eisenstein, priv. comm.). All magnitudes 
and colours presented t hroughout this pap er are corrected 
for Galactic extinction (|Schlegel et al.|[l998t ). 



3 MATCHING SAMPLES 

Different techniques were employed to select LRGs in the 
SDSS and the 2SLAQ survey, resulting in intrinsic differ- 
ences between the properties of the LRGs in each sample 
(Figure O. In particular, the magnitude dependent colour 
cut used in the SDSS selection results in only the very red- 
dest galaxies being included in the SDSS LRG sample at 
fainter magnitudes. Therefore, if we wish to make a mean- 
ingful comparison of the evolution of LRGs with redshift we 
must make additional colour and magnitude cuts to ensure 
that we exactly match the samples from the two surveys. 

Following Paper I, we assume that the evolution of the 
LRGs stellar populations can be approximated by simple 
passive ageing. We therefore use the same models as Pa- 
per I to generate K+e corrections which are used to correct 
the observed magnitudes of each sample to a common frame. 
Paper I demonstrated that these models do not perfectly de- 
scribe the colour evolution of the LRGs because of inadequa- 
cies in the stellar population synthesis models. To minimize 
the magnitude of these corrections, Paper I restricted their 
LRG samples to tight redshift ranges at approximately z = 
0.2 and z = 0.55 where the u,g and r filters approximately 
map onto the g, r and i filters respectively. These same red- 
shift cuts are again applied to the samples used herein. 

In this paper we take two approaches to matching the 
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Figure 1. The ^'^{g — i) versus Mo.2r colour magnitude relation 
for SDSS main galaxies with 0.15 < z < 0.21, all K+e corrected 
to 2; = 0.2. The black points in each panel show the whole sample. 
The top panel shows those galaxies that are selected to be in the 
2SLAQ selection matched sample (see text) when K-|-e corrected 
to z = 0.55 (green points). The bottom panel shows those galaxies 
that would be selected in the the SDSS selection matched sample 
(see text) when K+e corrected to z = 0.2 (red points). 



selection between these two redshifts. In the first we follow 
the procedure of Paper I. We take all the SDSS LRGs with 
0.17 < z < 0.24 and K+e correct their magnitudes to both 
z=0.2 and z=0.55. We then apply the SDSS selection crite- 
ria using the z = 0.2 magnitudes and the 2SLAQ selection 
criteria using the z — 0.55 magnitudes. We then execute the 
same procedure on the 2SLAQ LRGs within 0.5 < z < 0.6. 
We note that since the 2SLAQ selection is significantly bluer 
in the rest-frame than the SDSS selection; it is the applica- 
tion of SDSS selection cuts that is removing the majority of 
the LRGs removed from each sample by this procedure. We 
will therefore describe these samples as the SDSS selection 
matched samples. 

Our second approach makes use of the MAIN galaxy 
sample from the SDSS rather than just the LRG sample, 
although there is considerable intersection over the redshift 
range we are considering here. We limit the MAIN galax- 
ies to 0.15 < z < 0.21 and then apply our K+e corrections 
to correct to both z = 0.2 and z = 0.55. For the galaxies 
at z = 0.21 the Vpet = 17.77 magnitude limit of the MAIN 
galaxy sample corresponds to Mo. 2^ = -22.3. Mo. 2^ is cal- 
culated by determining the apparent magnitude the galaxy 
would have ai z — 0.2 in the SDSS r-band filter using our 
assumed K+e corrections, and is then converted to an abso- 
lute magnitude using the distance modulus without the use 
of any further K or evolutionary corrections. The Mo. 2^ — 
-22.3 is only 0.3 magnitudes brighter than the limit of the 
2SLAQ sample when K+e corrected to this redshift. Since 
the MAIN sample contains galaxies of all colours we can 
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Table 1. The redshift range, selection and number of galaxies in 
each sample defined in the text. 



Sample 


Redshift 


Selection 


Number 


1 


0.17 <z < 0.24 


SDSS 


9,912 


2 


0.5 < ^ < 0.6 


SDSS 


1,239 


3 


0.15 < z < 0.21 


2SLAQ 


11,350 


4 


0.5 < ^ < 0.6 


2SLAQ 


2,814 



generate a sample matching the 2SLAQ selection and we 
only need limit the 2SLAQ sample by this Mo. 2^ cut. There 
is, however, an additional complication. As shown in Paper 
I the errors on the photometry at the faint magnitudes of 
2SLAQ result in a large scatter of objects across the colour 
and magnitude selection boundaries. To mimic this effect 
we measure the magnitude error distributions of the 2SLAQ 
galaxies as a function of magnitude and modify the mag- 
nitudes of the SDSS galaxies randomly following this error 
distribution. We then apply the 2SLAQ selection criteria to 
both samples K+e corrected to z — 0.55 along with a cut at 
Mo. 2^ = -22.4. This slightly brighter cut than the Mo.2r = 
-22.3 limit allows the inclusion fainter galaxies which are be- 
ing scattered into the selection region by the application of 
the 2SLAQ photometric errors mimicking the effect present 
in the 2SLAQ data. We will refer to these samples as the 
2SLAQ selection matched samples. 

We are unable to account for the effect of the photomet- 
ric errors on the selection in the SDSS selected LRG sample 
as we don't have galaxies in that sample which are fainter 
or bluer than the LRGs. In Paper I we corrected the LP at 
z — 0.55 for this sample using a sub-region that had deeper 
photometry. We are unable to apply such a correction in this 
work since the significantly smaller area (~l/3 of the total) 
of this sub-region would result in a very poor measurement 
of the clustering and render any correction highly unreliable. 
The smaller area was not a problem for the LP measurement 
since the region of the LP most affected was the faint end 
were the galaxies were most numerous. The correction was 
also only required for a subsection of the LP, which one 
could always choose to disregard, whereas it would affect 
the entire correlation function. Por this reason, when mak- 
ing direct evolutionary comparisons between redshifts, we 
will focus on the 2SLAQ selection matched samples. Table 
[1] gives the number of galaxies in each sample and Pigure[T] 
illustrates the difference between the two selection criteria. 



4 THE MEASURED 2PT-CORRELATION 
FUNCTION 

The 2pt-correlation function, ^(r), is defined as the excess 
probability above Poisson of finding an object at a separa- 
tion r from another object. This is calculated by comparing 
the number of pairs as a function of scale in our galaxy cata- 
logues, with the number in a random catalogue, which covers 
the same volume as our data. W e make this measurement 
usme the lLandv fc Szalavl (|l993l ) estimator. 



1 

'RR 



DD 



-2Dr( — ) +RR 



(1) 



where DD, DR and RR are data-data, data-random and 
random-random pair counts respectively, and nn and na 
are number of galaxies in the data and random catalogues. 



4.1 Incompleteness corrections and error 
estimates 

When making the 2pt-correlation function measurement in 
our samples we must account for the varying completeness 
across our surveys. Por both the SDSS and 2SLAQ we sepa- 
rate the galaxies into unique regions based on the positions 
of the overlapping spectroscopic plates. Por the SDSS survey 
we use the regions defined in the SPSS C atalogue Archive 
Server fsee lAdelman-McCarthv et ai1l2007l . for details.). Por 
2SLAQ we use regions defined using the angular mask con- 
structed using repeated runs of the 2dP-configure software 
(see Paper I for details). Within each region we determine 
the number of targets with reliable redshifts [Nr) and the 
number of targets that could have been observed (Nt)- The 
completeness in each region is then defined as the ratio of 
these (Nr/Nt)- Regions with completeness below 65% are 
removed. 

To correct for the remaining incompleteness we wish 
to assign a weight ^ 1 to all the galaxies that have a re- 
liable redshift. We begin by assigning each target galaxy a 
weight equal to the inverse of the completeness of the re- 
gion in which it lies. Por those that do not have a redshift 
the weight is redistributed to its three nearest neighbours. 
This maintains some of the spatial information, although 
on the smallest scales where fibre collisions become impor- 
tant the clustering signal is likely to be underestimated. The 
weights of the galaxies with redshifts in a given region are 
then renormalised so that the mean weight in that region 
is as it was before the redistribution, i.e. the inverse of the 
completeness. 

An alternative approach, often used to correct for in- 
completeness, is to reduce the number of random points in 
regions with low completeness. We do not do this for two 
reasons. Pirstly, by having regions with lower numbers of 
random points we will be unnecessarily increasing the noise 
in these regions. Secondly, and more importantly, unlike in 
the SDSS, the spectroscopic plates in 2SLAQ were evenly 
spaced with no allowance made for the variation of the tar- 
get space density. This means that regions with a high target 
density (i.e. highly clustered regions) will be more likely to 
have a lower completeness. We calculate the completeness in 
regions defined by the overlapping plates and so by simply 
reducing the number of random points based on this com- 
pleteness we would be likely to systematically underestimate 
the clustering on scales smaller than the given region. We 
would be preferentially removing the most clustered galeixies 
and then renormalising the clustering calculated from the re- 
maining less clustered galaxies by the ratio of the number re- 
moved (i.e. the completeness). Since we instead redistribute 
the weight of the galaxies without redshifts to their nearest 
neighbours, we are likely to be up weighting other galaxies 
in the most clustered regions and will therefore be making 
a better estimate of the true clustering amplitude. 

Nearly all of the completeness regions have annular 
scales up to 2 degrees which corresponds to 32.6 Mpc 
at z — 0.55 and so this effect is likely to be important over 
nearly all the scales we consider in this paper. In fact the 
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Figure 2. The redshift space 2pt-correlation functions at z ^ 0.2 (red open circles) and z ~ 0.55 (blue filled circles) and their ratio [z < 
0.55/2 ~ 0.2) for the SDSS selection matched (left) and for the 2SLAQ selection matched (right) samples. 




Figure 3. The projected 2pt-correlation functions at 2 ~ 0.2 (red open circles) and 2 ~ 0.55 (blue filled circles) and their ratio (2 ~ 
0.55/2 ~ 0.2) for the SDSS selection matched (left) and for the 2SLAQ selection matched (right) samples. The lines show power law fits 
on scales 0.32 < rp < 32h~^ Mpc. 



clustering is ~ 5% lower for the 2SLAQ samples when cal- 
culated by just reducing the number of randoms. 

We generate random catalogues for each galaxy sample 
following the angular masks of the surveys with constant 
space density and 20 times the number of random points 
as data. The regions around bright stars are removed from 
both data and random catalogues, as galaxies in these re- 
gions are known to have systematically incorrect magnitudes 
due to poor sky subtraction in SDSS photometric pipeline 



('Mandelbaum et al.l 120051 : lAdelman-McCarthv et al.ll2006l ). 

Rcdshifts are assigned to the random catalogues by ran- 
domly sampling a polynomial fit to the redshift distribution 
of each galaxy sample. We note that within the tight red- 
shift ranges of the samples considered here all the samples 
are approximately volume limited. 

We estimate the errors on our 2pt-corr elation function 
measurements using jack knife re-sampling (jScranton et al.l 
l2002l : IZehavi et aLli2005h . We split the SDSS area into 40 
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equal area regions and the 2SLAQ area into 32 equal area 
regions. We then calculate each 2pt function removing one 
area at a time to generate a full covariance matrix. Through- 
out this analysis we measure the pair counts using the KD- 
tree c ode in the NTROPY software package l|Gardner et al.l 
I2OO7I ). 

4.2 Various clustering estimators 

The peculiar velocities of galaxies generate errors in the dis- 
tance measurements along the line of sight. This means that 
our basic measurement of ^, which is based on redshift dis- 
tances, is affected by these redshift space distortions. By 
separating the clustering signal into contributions perpen- 
dicular (rp) and parallel (tt) to the line-of-sight {S.{rp,n)) 
and then integrating over the tt direction, one obtains the 
projected correlation function 

w,ir,) = 2^°° dnarr,n) = 2^"" j/'^l, - (2) 

The final expression only involves the real space corre- 
lation function ^(r) showing that Wp(r„) is not compro - 
mised by redshift space distortions l|Davis fc Peeble3ll983l ). 
One can invert equation ^ by interpolating between the 
binned w{rp) to yield an e stimate of ^(r) whic h is free of 
redshift space distortions ([Saunders et al.lll992l '). If ^(r) = 

i r/rp) then equati on Q can be solved analytically 
Davis fc Peeblesll 19831 ). 

In practise, one models Wp with the second of the equal- 
ities above, but measures it using the first. However, when 
making the measurement, it is only sensible to integrate out 
to some maximum n because ^{rp,n) is poorly known on 
very large scales. We integrate to 80/i~^Mpc which appears 
to give stable results. 



4.3 Observed evolution of clustering 

Figures [2H11 show ^(s), and ^(r) for the four samples 

described in Section [31 along with the ratio of the functions 
between the two redshifts. Figures [SHU also show the result 
of fitting power laws over the scales 0.32 < rp < 32h~^ Mpc 
using the full covariance matrices derived from the jackknife 
re-sampling technique. We limit the fits to scales greater 
than 0.32/i~^ Mpc since our weighting scheme does not fully 
correct for the effect of fibre collisions on smaller scales. Ta- 
ble [2] provides the best fit values of ro, 7 and the associated 
reduced x^j with error contours shown in figure [S] These 
measurements show that there is very little evolution in the 
clustering amplitude of LRGs between z ~ 0.55 and z ^ 0.2, 
but there is a marginally significant increase in the slope. 



4.4 Comparison with previous work 

Several previous studies have performed similar analyses to 
those we present here; it is important to make a compari- 
son of the results befo re further investigat ing the meaning 
of these measurements. IZehavi et al.l (|200^ present the 2pt- 
correlation function for three slightly different samples of 
SDSS LRGs. One of these samples, with -23.2 < Mg < 
—21.2, has an almost identical space density to the z; ~ 0.2 



SDSS selection matched sample, although at a higher red- 
shift = 0.28). The 2pt-correlation functions in redshift, 
projected and real space for this sample are almost indistin- 
guis hable, with i n the errors, to those presented here. 

IRoss et al.l (|2007l ) present measurements of the 2pt- 
correlation function for the Sample 8 2SLAQ LRGs. This 
sample is similar to the 0.5 < z < 0.6 2SLAQ selec- 
tion matched sample although with a larger redshift range 
and slightly fai nter absolut e mag nitude cut. The power-law 
fit to w{rp) in IRoss et all (|2007i ) has a very similar slope 
(7 — 1.83 ±0.05) to that measured here with a lower ampli- 
tude (ro = 7.3 ± 0.34fe~^ Mpc) . This lower amplitude is to 
be expected as I Ross et al.l (|2007t ) include intrinsically fainter 
galaxies in their sample. To make a more direct comparison 
we recalculat ed w(rp) us i ng a s election almost identical to 
that used bv IRoss et al.l (|2007l ). This produces an almost 
identical slope (7 = 1.81 ±0.03) but a slightly hi gher ampli- 
tude (ro = 7.85±0.15/i"^ Mpc) to that found bv lRoss et all 
(|2007l ). This is to be expected, as IRoss et "all (|2007^ sim- 
ply reduced the number of randoms points as a function of 
completeness. As discussed in Section |4l the 2SLAQ data 
are more likely to be incomplete in the densest regions and 
so by reducing the number density of random points as a 
function of completeness they will tend to underestimate 
the clustering on scales smaller than the regions in which 
the determine the completeness. 

Since we produ ce almost id e ntical measu rements to 
those presented in IZehavi et al.l l|2005h and IRoss et"aLl 
((20071 ) ■ with a completely independent analysis and differ- 
ent techniques on largely the same data, we can be confident 
that our measurements are accurate. We now consider what 
our measurements imply for our LRG samples. 

The slope of the 2pt-function is known to depend on 
colour/sp ectral type: bluer galaxies have a s hallower slope 
(e.g. Zch avi et al.l l2002l : iNorberg et al.ll2002l ). Could it be 
that there are more blue galaxies in the z — 0.55 samples? 
We have assumed passive evolution when defining the sam- 
ple selection, so it seems unlikely that this would include 
more intrinsically bluer/later- type galaxies at high redshift 
than at low redshift. We could, however, be scattering more 
blue galaxies across the selection boundaries at high red- 
shift than at low redshift, for instance, if there were more 
galaxies populating the blue cloud close to the red sequence 
at high z. If this is the case one might expect to see a dif- 
ference in the slopes between the 2SLAQ selection matched 
and SDSS selection matched samples, as the SDSS selec- 
tion only allows the reddest galaxies to be included at the 
faintest magnitudes where the scattering is most significant. 
This is not the case, suggesting that despite the fact that we 
have selected galaxy populations consistent with purely pas- 
sive evolution, both dynamically and in terms of their stellar 
populations, we are in fact seeing some additional evolution 
in the LRG population. 



4.5 Comparison with a no-merger model 

If, as suggested in Paper I, the LRGs do not merge with one 
another, then the large scale bias is predicted to evolve as 
bio = 1 ± (bhi — i){phi/Dio) where D is the linear growth 
factor (|Mo fc Whitell 19961 : [jvvll 19961 ). In this case, the ratio 
of the correlation functions should be 
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Figure 4. The real space 2pt-correlation functions at z ~ 0.2 (red open circles) and z ~ 0.55 (blue filled circles) and their ratio (z 
0.55/2 ~ 0.2) for the SDSS selection matched (left) and for the 2SLAQ selection matched (right) samples. The lines show power law fits 
on scales 0.32 < r < 32h-^ Mpc. 




Figure 5. 68%, 90% and 99% confidence intervals for power law fits on scales 0.32 < rp < 32h^^ Mpc to the projected 2pt-correlation 
functions (left) and the real space 2pt-correlation functions (right) at z ~ 0.2 (red) and z ~ 0.55 (blue) for the SDSS selection matched 
(dashed lines) and for the 2SLAQ selection matched (solid lines) samples. The error bars show the 1 a errors on the individual parameters. 



C<o(r) blDl y ho ) 

on large scales. Note that this diflers from the growth of 
the dark matter clustering strength, because of the factor 
(bhi/bio)^. Since Dm/Dio ^ 1, the large scale clustering 
strength should increase at late times. For Zhi ~ 0.55 and 
zio — 0.2 in our chosen cosmology, Dhi/Dio — 0.84. We will 



argue below that bhi/bio — 2.16/1.91 — 1.13, so the expected 
ratio of large scale clustering strengths is 0.9. 

A similar argument can be ma de for the clu stering in 
redshift-space: on scales where the iKaised l| 19871 ) analysis 
of redshift space distortions applies, the expected ratio of 
redshift-space clustering amplitudes is 

_ l + 2/3,„/3 + /3g,/5 / b^OM X' 
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Table 2. Values of the power-law fits and the reduced to w(rp), and ^(t) in the range 0.32 < r < 32h-^ Mpc. 



Selection 


Redshift 


ro (h 


1 Mpc) 


7 




- 






w{rp) 




w{rp) 




w{rp) 




SDSS 


0.21 


9.47 ± 0.29 


9.52 ± 0.39 


1.96 ± 0.03 


1.87 ± 0.04 


0.76 


0.54 


SDSS 


0.55 


9.61 ± 0.62 


9.42 ± 0.76 


1.79 ± 0.06 


1.73 ± 0.09 


0.53 


0.71 


2SLAQ 


0.19 


7.64 ± 0.29 


7.72 ± 0.36 


1.98 ± 0.04 


1.89 ± 0.05 


1.58 


1.10 


2SLAQ 


0.55 


8.29 ± 0.30 


8.15 ± 0.42 


1.77 it 0.06 


1.71 ± 0.08 


1.02 


0.80 



where Pio ~ nf/^/&;o and (3^ ^ fil^^* /bm- Again, in the no 
merger model, the low redshift population is expected to be 
more strongly clustered. For the two LRG samples studied in 
the main text, the expected ratio is (1.236/1.211)0.9 — 0.92. 

Figures [2J [3] and [3] show the ratios of ^(s), Wp and 
^(r) measured at z ^ 0.55 and z ~ 0.2 which appear to 
be consistent with little or no evolution. The expected ra- 
tios calculated above are inconsistent with the data at the 
93% level for ^(s) and the 80% level for Wp on large scales 
(r > 3h~^ Mpc) where these calculations apply. Thus, the 
clustering signals suggest that the low redshift LRG popu- 
lations are not simply passively evolved versions of the high 
redshift population, although we are not able to conclusively 
demonstrate this with the large-scale clustering measure- 
ments alone. In the following sections we model both the 
evolution of the clustering on all scales and the number den- 
sity to further constrain the evolution of LRGs. 



5 HALO MODEL ANALYSIS 

The halo model (see Cooray & Sheth 2002 for a review) as- 
sumes that the galaxy clustering signal encodes information 
about the Halo Occupation Distribution (HOD) - how the 
galaxies populate Dark Matter haloes - in particular, how 
the HOD depends on halo mass. This approach has recently 
been used to constrain the HODs of galaxies in a number of 
large data sets. We apply such a model here to try to gain 
insight into our LRG populations, how they have evolved, 
and how well or otherwise this evolution can be described by 
the passive no-merger model. Our analysis of the no-merger 
model has strong s imil arities to that rec ently performed by 
IWhite et af] (|2007l ) and lSeo et all (|2007l ). However, whereas 
their work was primarily numerical, our analysis shows that 
the entire discussion can be analytic. 



5.1 The centre-satellite HOD 

In the halo model, every galaxy is associated with a halo; all 
haloes are 200 times the background density whatever the 
mass M of the halo. Sufficiently massive haloes typically 
host more than one galaxy. The halo model we use distin- 
guishes between the central galaxy in a halo, and the oth- 
ers, which are u sually called satellites . This is motivated by 
simulations (e.g. iKravtsov et al.ll200j ). and has been a stan- 
dard assumption of ^emi-analytic galaxy formation models 
for many years (e.g. iBaugh.2006. ) . There is now strong obser- 
vational evidence that the two types of galaxies are indeed 
rather different, and that the halo model parametrisation of 
this d ifference is rather accurate (jSkibba. Sheth fc Martiriol 
l2007h . 



The fraction of haloes of mass M which host centrals is 
modelled as 



{Nc\M) = exp{-Mrr^^„/M). 



(5) 



Only haloes which host a central may host satellites. In such 
haloes, the number of satellites is drawn from a Poisson 
distribution with mean 

{Ns\M} ^{M/Khr. (6) 
Thus, the mean number of galaxies in haloes of mass Al is 
{N\M) = {iVc|M)[l + {N,\M)], (7) 
and the predicted number density of galaxies is 



Ug = / dMn{M) {N\M), 



(8) 



where n{M) is the halo mass function, for which we use the 
parametrisation given by Sheth & Tormen (1999). 

We further assume t hat the satellite galaxies in a 
halo trace an NFW profile (|Navarro. Frenk. fc Whitelll996l ) 
around the halo centre, and that the haloes are biased trac- 
ers of the dark matter distribution. The halo bias depends 
on halo mass in a way that can be estimated directly from 
the halo mass function (Sheth & Tormen 1999). With these 
assumptions the halo model for ^(r) is completely specified 
(e.g. Cooray & Sheth 2002). We then calculate w{rp) from 
^ using the second of equations ([2]). 

In addition to 5, we are interested in the satellite frac- 
tion. 



Fsat = / dMn{M) (N,\M) {Ns\M)/ng 



(9) 



and two measures of the typical masses of LRG host haloes: 
an effective halo mass 



Me// = J dMMn{M) {N\M)/ng, 
and the average linear bias factor 



bg= I dMn{M)b{M){N\M)/ng, 



(10) 



(11) 



where b{M) is the halo bias. 

Our notation is intended to make explicit the fact 
that the mean number density of central-satellite pairs 
from such haloes is n{M) {Nc\M) {Ns\M) , and the mean 
number density of distinct satellite-satellite pairs is 
n{M) {Nc\M) {Ns\Mf /2 (because we are assuming the 
satellite counts are Poisson). 

For completeness, our model for the real-space 2-point 
function is 



(12) 
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Figure 6. HOD fits on scales 0.32 < rp < 50h ^ Mpc to the projected 2pt-correlation functions at z ^ 0.2 (red) and z ^ 0.55 (blue) for 
the SDSS selection matched (left) and for the 2SLAQ selection matched (right) samples. 



where 



5.2 HOD fits 



- /dM!^Wi^(7V.|M)^ (13) 



and 

with 

P2h{k) 



n(M){N,\M) {N^Mf \{r\M) 

Hg 2 ngM2 ' 



f dkk^ 
J k 



= bg{kf PLin{k), 



where 



(15) 
(16) 

1 + {Ns\M)u{k\M) 



In the expressions above, p{r\AI) is the density profile of 
haloes of mass M, X{r\M) denotes the convolution of two 
such profiles, u{k\M) is the Fourier transform of p{r\M)/M , 
and Phin{k) denotes the linear theory power spectrum. In 
practise, we usually approximate bg{k) by its value bg aX k = 
(equation llip . All these quantities, along with the mass 
function n{M) and bias factor b{M), are to be evaluated 
at the redshift of interest. We have already specified how, 
for a given halo mass, the virial radius depends on redshift; 
the NFW halo density profile is also specified by its con- 
centration, which w e assume is c = 9 (M/Mto)~°'^^/(l + z) 
l|Bullock et aLllioOll ). All this, in the right hand side of equa- 
tion ([2|, gives the halo model calculation of Wp{rp). 

Our halo model calculation of ^(s) makes two addi- 
tional assumptions: first, that satellite galaxies within haloes 
have isotropic velocity dispersions which are proportional to 
GM/rvir, and second, that the motion of the centre of mass 
of a halo is well described by linear theory. 



We fit for the parameters Mmin, Mi and a (see equations [5] 
andE]) by minimising a defined as the sum of the squared 
difference between the predicted and measured Ug and w(r'p) 
for a range of rp. We use w{rp) rather than ^(r) as the 
numerical inversion required to calculate ^(r) increases the 
uncertainties and systematically reduces the slope in our 
power law fits. Our fitting makes use of the full covariance 
matrices over 0.32 < Vp < 50 Mpc. We exclude scales 
smaller than 0.32 Mpc as we are not confident that we 
have sufficiently corrected for fibre collisions. We note that 
the best fitting parameters are not significantly changed if 
the smallest bin included in the fit is one smaller or larger. 

The errors on the fits are determined by finding the re- 
gion of parameter space with a Sx^ ^ 1 (la for 1 degree of 
freedom) from the best fit and then determining the max- 
imum and minimum parameter values within that region. 
For biin, Meff, and Fsat, which depend on all three of the 
fit parameters, the region used contains Sx'^ ^ 3.53 (la for 
3 degrees of freedom). 

The resulting best fits are shown in Figure [6] and the 
best fit values for the HOD parameters are given in Table |31 
We have checked that our best fitting model also provides 
a good description of our measurements of ^(s) and S,{r). 
These parameters were not included in our definition of x^ 
because the halo model of ^(s) requires further assumptions 
than does w{rp). Table|3]also provides the associated values 
of Fsat, Mcff and fcun. 

The best fitting HODs are shown in Figure[7] Increasing 
ag (see Table I?} increases Mmin and Mi, and decreases a. 
The bias decreases to compensate for the increased cluster- 
ing strength of the dark matter, and Me/ / increases because 
Mmin is larger. The satellite fraction remains approximately 
the same, as a has reduced to compensate for the increase 
in Ml. 
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Table 3. The best fitting HODs to Wp(rp) assuming eg = 0.8 

Selection Redshift Density M^in Mi a xled ^eff Fsat 

(IQI^Mq) (IO^A/q) (IQI^Mq) {%) 

SDSS 0.21 0.94 ± 0.01 3.80 ± 0.07 34.2 ± 2.1 1.67 ± 0.23 1.22 2.11 ± 0.03 9.52 ± 0.59 10.1 ± 3.7 

SDSS 0.55 0.73 ± 0.02 3.46 ± 0.06 34.0 ± 2.5 2.10 ± 0.38 1.12 2.42 ± 0.05 6.24 ± 0.51 4.7 ± 2.5 

2SLAQ 0.19 1.64 ± 0.01 2.44 ± 0.02 27.0 ± 1.1 1.58 ± 0.13 0.77 1.91 ± 0.02 7.62 ± 0.41 10.4 ± 2.1 

2SLAQ 0.55 1.65 ± 0.03 1.88 ± 0.02 21.8 ± 1.5 2.02 ± 0.2 1.23 2.16 ± 0.03 4.76 ± 0.20 6.2 ± 2.3 



Table 4. The best fitting HODs to Wp(rp) assuming eg = 0.9 

Selection Redshift Density Mmin Mi a x^ed ^eff Fsat 

(10-4?t3Mpc-3) (IO^^Mq) (IO^Mq) (IO^^Mq) (%) 

SDSS 0.21 0.94 ± 0.01 4.43 ± 0.15 45.5 ± 5.7 1.38 ± 0.16 1.45 1.91 ± 0.03 11.82 ± 0.60 11.8 ± 2.2 

SDSS 0.55 0.73 ± 0.02 4.15 ± 0.09 46.3 ± 3.9 1.91 ± 0.39 1.13 2.20 ± 0.04 8.22 ± 0.81 5.1 ± 2.8 

2SLAQ 0.19 1.64 ± 0.01 2.77 ± 0.03 34.2 ± 1.3 1.38 ± 0.13 0.97 1.73 ± 0.02 9.59 ± 0.68 11.7 ± 2.3 

2SLAQ 0.55 1.65 ± 0.03 2.19 ± 0.03 28.2 ± 2.1 1.86 ± 0.20 1.26 1.96 ± 0.02 6.28 ± 0.36 6.8 ± 2.4 



For our standard choice of as = 0.8, the LRGs populate 
haloes with masses of order lO" - IO^Mq; most of these 
LRGs are central galaxies - the satellite fractions are typi- 
cally less than 10%. In the lower redshift samples Meff is 
larger by about 50%, the bias is smaller by about 10%, and 
the satellite fraction has approximately doubled. The growth 
in Msf / is a consequence of a 10% increase in Mmin, a small 
decrease in Mi /Mmin, and a significant decrease in a. 

It might seem paradoxical that decreasing a increases 
the satellite fraction. This is a consequence of the fact that 
Ml is larger than the mass scale on which the halo mass 
function drops exponentially (for as = 0.8, this scale is 
0.6 X W^^h-^MQ and 1.9 x W^^h'^MQ a,t z = 0.55 and 
z — 0.2 respectively; when erg — 0.9, these masses become 
1.3xl0^^h-^M@ and 3.9 x IO^^/i'^Mq). Thus, increasing a 
increases the number of satellites in (the exponentially rare) 
haloes more massive than Mi but decreases the number in 
less massive haloes which are exponentially more abundant. 

The larger satellite fractions at low redshift are best 
understood by thinking of the central and satellite popu- 
lations separately. If there is no merging, then the high z 
satellites are satellites even at low z, whereas some of the 
high z centrals have become satellites at low z (e.g., if their 
host halo merged with a more massive halo). As a result, the 
satellite fraction increases. Merging would act in the oppo- 
site sense (satellites merging with satellites or with centrals 
would both reduce the satellite fraction) . 

We note that the best fitting HODs for the z = 0.55 
samples are in excellent agreeme nt with those presented in 
iBlake. CoUister. fc Lalia^ (|2008f ) who fit HODs to the an- 
gular clustering of 380,000 LRGs selected using the 2SLAQ 
LRG selection criteria with photometric redshifts 0.45 < 
Zphot < 0.65. 



6 CONSTRAINING LRG MERGERS 

Paper I demonstrated that the evolution of the LF of LRGs 
was consistent with passive evolution of the stellar popu- 
lations, and did not require any merging. If true, then as 



discussed in Section 14.51 the bias should evolve as b{zio) = 
l + (b(zhi) — l)p(zhi) / D(zio) w here D(z) is the growth factor 
(|Mo fc Whitll 19961 : IfVvII 19961 ). When applied to the bias of 
the best fitting z = 0.55 HODs for the two samples, the pre- 
dicted bias factors are 1.98 ± 0.02 at z ^ 0.19 for the 2SLAQ 
selected sample and 2.20 ± 0.04 a.t z = 0.21 for the SDSS 
selected sample. Both these values are significantly larger 
than the measured values given in Table O with the evolu- 
tion in the 2SLAQ selected sample bias being incompatible 
with no-merging hypothesis at a significance of 98.4%. This 
is at a higher significance level to that calculated in Sec- 
tion 14.51 using just the ratio of the large-scale clustering; 
the inclusion of the number density constraints in the HOD 
fits results in significantly smaller relative errors on the bias 
measurements than would be derived using clustering alone. 

This argument against pure passive evolution still uses 
only the large-scale clustering signals at the two epochs. 
In what follows, we use the language of the halo model to 
show that the evolution of the small-scale clustering signal 
also contains interesting information, and can provide even 
greater constraints on the importance of merging. 



6.1 HOD Evolution: No mergers 

If we specify how galaxies populate haloes at some early 
time, (TVjm), then we can estimate how this evolves as the 
haloes merge. If the haloes merge but the galaxies do not, 
then 

(N\M) = / dmN{m\M) (TVjm) = C(M) + S{M) (17) 
Jo 

where N{m\M) is the mean number of haloes of mass m 
which are in haloes of mass M at the later time, and 

C(M) = / AmN{m\M){Nc\m) and (18) 

I'M 

S{M) = / AmN{rn\M){Nc\m) {N,\m); (19) 

^0 
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Figure 7. The mean number of LRGs per halo as a function of halo mass (top) and the mean number of LRGs per halo times the 
number density of haloes as a function of mass (bottom) at ^ ~ 0.2 (red) and z ~ 0.55 (blue) for the SDSS selection matched (left) 
and for the 2SLAQ selection matched (right) samples. The total, central and satellite contributions are shown by the solid, dashed and 
dotted lines respectively. 



For N{m\M) we use the expressions given by Sheth & Tor- 
men (2002), which generalise those of Lacey & Cole (1993). 
The Appendix shows that this guarantees that the comoving 
density Ug is constant, whereas the large scale bias evolves 
in accordance with the continuity equation. 

Whereas C{M) counts the objects which used to be 
centrals, S{M) counts the satellites. Note that although 
(Nc\m) ^ 1, there is no guarantee that C{M) ^ 1; indeed, 
for M Mmin, one expects C{M) ^ L Figure [S] shows 
this explicitly; at late times, massive haloes may host many 
galaxies which were centrals at the earlier time. 

If we force {N\M} to have the same functional form as 
{N\m}, then we can fit for Mmin, Mi and a at the later 
time. These fitted values can then be inserted into the halo 
model calculation of ^. The Figure |8] shows that forcing this 
parametrisation allows a good but not perfect description of 
the passively evolved HOD: the passively evolved HOD has 
a more gradual transition from to 1. 

It will turn out that, for the present study, it is impor- 
tant to accurately model this transition. This is because we 
are studying rare objects which populate the high-mass end 
of the mass function. As a result, haloes which host zero 
or one galaxies are substantially more numerous than those 
which host more. Hence, allowing some lower mass haloes to 
host more than one galaxy (while making more such haloes 
void of galaxies) can affect the number of small separation 
pairs substantially. 

To illustrate this effect, let po{M) denote the probabil- 
ity that a halo of mass M contains no galaxies which were 



centrals at the higher redshift. Then 

(iVe|A/> = l-po{M) and (20) 
{N,\M){Ns\M} = S{M) + C{M)- {Nc\M). (21) 

The second equation assumes that only one of the high-z 
centrals in a halo continues to count as the low-z central; 
the others (of which there are C{M) — {Nc\M) on average) 
count as \ow-z satellites. The mean galaxy count {N\M) 
is given by inserting these expressions in equation ((Tjl. This 
exercise shows that the problem is to model po{M); the next 
subsection studies three different models. 



6.2 HOD Evolution: Small scale clustering and 
the abundance of empty haloes 

The quantity po{M) counts the number of haloes of mass 
Af which were formed from mergers of objects which con- 
tained no galaxies. If the threshold A/min were sharp, then 
this would be simply related to the number of haloes at low 
redshift which did not have a single high-redshi f t prog enitor 
of mass greater than Afmin. ISheth fc Lemsor] (|l999l ) have 
studied this problem; they provide expressions for the fc-th 
factorial mom ent fj,k of t he progenitor distribution. (Results 
in [ciis as- Miranda et al.l (2002) suggest that these expres- 
sions are quite accurate.) In principle, these can be used to 
estimate po, since po = 1 + i^k/k\, where the sum 

runs from fc = 1 to an upper limit which is set by mass 
conservation; a halo of mass M can have at most A//A/min 
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Figure 8. The mean number of LRGs as a function of halo 
mass ai z ^ 0.19 obtained by passively evolving the best fit- 
ting z = 0.55 2SLAQ HOD to ^ = 0.19. Upper dashed and lower 
dotted lines show the contributions from objects which used to 
be centrals and satellites; they sum to give the solid curve which 
drops to zero at smaller mass scales; lower dashed and upper dot- 
ted lines, which sum to give the other solid curve, show the result 
of fitting this (N\1\I) to the form given in equations JSj and (|6j. 



progenitors. If the HOD where a step function then Mmin 
would be the same as in Equation [5] else, it need not be. In 
practise, this is a complicated sum, so we have studied a few 
simpler models. 

In our first model, we set 

-C{M) 



Po{M) 



(22) 



This would be appropriate if the distribution of the number 
of high-redshift centrals in low-redshift haloes were Poisson 
(so /ife — fii), with mean fii = C{M), and only one of these 
centrals continues to count as the low-redshift central; the 
others count as low-z satellites. Note that if C <^ 1, then 
{Nc\M) — > C(M), so there is no correction to the satellite 
counts. And if C > 1 then (A^c|M) 1 and the satellite 
counts are increased by C — 1. Thus, our model interpolates 
smoothly between these two sensible limits. We show the 
resulting evolution in the way galaxies populate haloes and 
in the clustering in the top panels of Figures |9] and [10] as the 
red lines. 

We have also studied what happens if, instead, we 
require a sharp transition between these two limits: set 
{N^\M) = C{M) and {Ns\M) = S{M) when C(M) 1, 
and {Nc\M) = 1 and {Ns\M) = S(Af)-hC(M)-l otherwise. 
Compared to the Poisson model, this model has many more 
low-redshift haloes which host a single central high-redshift 
galaxy, and few which host more than one such galaxy; the 
Poisson model has fewer haloes which host galaxies, each 
allowed to host more than one high-redshift central. This 
decreases the number of high-redshift central pairs in haloes 
(compared to the Poisson model), which means that the 
number of central-satellite pairs is decreased, thus decreas- 



ing the small scale clustering signal. (Of course, higher-order 
statistics will also be affected: the probability of finding a 
large region devoid of galaxies will be larger in the Poisson 
model.) This model is plotted as the blue lines in the top 
panels of Figures [9] and 1101 

Whereas this second model is perhaps too simple, the 
Poisson model almost certainly allows too many low mass 
haloes to contain more than one galaxy, thus resulting in 
too ma ny small scale pairs. Indeed, mass conservation argu - 
ments (jSheth fc LemsonI 1 19991 : ICasas-Miranda et all 12002( 1 
strongly suggest that the progenitor counts should be sub- 
Poisson (/ife < pi), especially at low masses. Furthermore, 
sub-Poisson counts are clear ly seen in the numerical models 
10 and 30 of lSeo et al.l l|2007f ). The following Binomial model 
conserves mass, and lies between these two extremes: 



Po(M) 



C{M) 



(23) 



where A'^max = int(Af/Mniin). We use this model as written 
for illustrative purposes only: in reality Mmin is unlikely to 
be the same quantity as in equation ([Sjl, and the integer 
changes in A^max as M increases produce artificial disconti- 
nuities in {N\M). Nevertheless, this model predicts a small 
scale clustering signal which lies below that associated with 
the Poisson model, but above that for the sharp threshold 
model shown as the green lines in the top panels of Figures 
Handlini 

The top panel of Figure[TD]shows that for all three mod- 
els for po{M) the passive evolution of the clustering predicts 
a far greater increase in the clustering strength than is ob- 
served. This is caused by the presence of too many satel- 
lite galaxies, with satellite fractions of 27±3%, 11±1% and 
19±2% for the Poisson, Step and Sub-Poisson models re- 
spectively compared to 10±2% for the best fitting HOD to 
the data. 



6.3 HOD Evolution: Central-central mergers 

Once we have decided how likely it is that a low-redshift halo 
contains at least one high-redshift central galaxy, we also 
study models in which centrals merge onto centrals. This is 
motivated by the fact that central galaxies are expected to 
be more massive than satellites, so dynamical friction may 
be more effective at making these objects merge onto the 
true low-redshift central. To model this case, we again use 
equation (1201) for {Nc\M), but we set 



{N,\M){Ns\M) = S{M)+f, 



no — mcrec 



C{M)-{Nc\M) ,(24) 



where /no-mcrgo is the fraction of low-redshift satellites 
which were high-redshift centrals, and have not merged with 
one another or onto the new central object. 



When /, 



no — merge 



1 then this is the same as the no 
merger model of the previous section; when /no-morgo = 
0, then the central galaxies of all the high-redshift haloes 
which merged to make a low-redshift halo have merged to 
make a single massive central galaxy. Strictly speaking, the 
model says nothing about what these objects merged with 
- they may have merged with one another or with other 
satellites - it only assumes that the number of objects which 
merge scales with M in the manner given above. However, 
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Figure 9. The mean number of LRGs per halo as a function of 
halo mass at 2 = 0.19 for the 2SLAQ selection matched samples. 
The top panel shows the effect of passively evolving the z = 0.55 
HOD to z = 0.19 using the three models for po{M) along with 
the measured HOD from the z = 0.19. The effect of including 
merging of the central galaxies for the same models is shown in 
the bottom panel. The total, central and satellite contributions 
are shown by the solid, dashed and dotted lines, respectively. 



the assumption that they merged onto the central object has 
considerable physical appeal. 

The result of applying this merger model for the three 
parametrisation of pq{M) are shown in the bottom panels 
of Figures [9] and 1101 For each model we chose the value of 
fno-merge that bcst matchcs the large scale clustering, 0.1, 
0, 0.25 for the Poisson, Step and Sub-Poisson models respec- 
tively. In all cases the agreement in the high mass haloes is 
much improved and the satellite fraction reduces to 11 ±3%, 
7±2%, 10±3%, comparable to the meaured value. The best 
fit at small scales is provided by the sub-Poisson model; this 
is reasuring, as it is the most physically motivated - although 
our implementation is not yet ideal. This suggests that the 
data are consistent with a generic prediction of hierarchical 
models - that the scatter in merger histories should pro- 
duce sub-Poisson scatter. The step model produces far too 
little small scale clustering, consistent with its lower satel- 
lite fraction, with both the Poisson and sub-Poisson models 
providing a reasonable match within the errors. 

We show in Figure[TT]a more detailed comparison of the 
passive and merger sub-Poisson model with the measured 
correlation functions by dividing each by the best fit to the 
z = 0.19 measurement. Also shown are the la confidence 
regions calculated by propagating the error on the fit at z 
— 0.55. This figure explicitly shows that the passive model 



o 




tp (h 'Mpc) 

Figure 10. The projected correlation function resulting from the 
evolving the 2 = 0.55 2SLAQ selection matched sample HOD 
to 2 = 0.19 using the three models for po(M) along with the 
correlation function from the HOD fit to the 2 = 0.19 sample. 
The the effect of including merging of the central galaxies for the 
same models is shown in the bottom panel. 



is ruled out at high significance. On large scales (> 3 h~^) 
Mpc the passive model is incompatible with the measured 
clustering at 2 = 0.2 at the 98% level, consistent with the 
constraints from the bias evolution given above. However, 
when smaller scales are included the passive model becomes 
increasingly incompatible with the measured clustering; for 
scales larger than 1 Mpc the passive model is excluded 
at a confidence level of greater than 99.9%, with the level 
of significance increasing with the inclusion of even smaller 
scales. The sub-Poisson merger model is consistent with the 
data on all scales, even though the fraction of centrals which 
are alowed to merge is determined by matching only the 
large-scale clustering. 

We have demonstrated that it is necessary to allow some 
merging (or some other method of removal) of some fraction 
of the high redshift LRGs if we wish to reproduce the clus- 
tering at low redshift. This will have the effect of reducing 
the space density of the evolved population at low-redshift, 
something that we do not observe in the data. The change 
in the space density associated with the best fitting sub- 
Poisson model is 9.2±2.6%, suggesting that at most about 
20% of the LRGs are merging with each other. In fact there 
are on average 2.34 high-redshift centrals in each merged 
halo, resulting in 16.1±4.6% of LRGs experiencing an LRG- 
LRG merger. This is consistent with the constrainst pro- 
vided by the luminosity function evolution of Paper I. For 
comparison the Possion model predicts a change in the space 
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Figure 11. Halo model fits to w{rp) at z ~ 0.19 (red solid) and 
2 ~ 0.55 (blue solid) for the 2SLAQ selection matciied sample. 
The effect of passively evolving the z = 0.55 HOD to z = 0.19 is 
shown as the dotted line and the effect of including merging of the 
central galaxies is shown as the dashed line. The bottom panel 
shows the ratios of the w(rp) fits shown above to the measured 
z = 0.19 fit. The shaded areas enclose the 1 sigma confidence 
regions. 



density of 19.4±5.5%, suggesting that up to 40% of LRGs 
have been involved in a LRG-LRG merger. This number is 
highly inconsistent with the LF measurments and lends fur- 
ther support to the sub-Poisson model. 

If we continue with the hypothesis that LRGs are merg- 
ing with one another, it is reasonable to assume that some 
red galaxies too faint to be included in our sample at z — 
0.55 will have also merged by 2 = 0.19, some of which will 
now be sufficiently luminous to be included in that sam- 
ple. These galaxies will then increase the space density of 
the low-redshift LRG sample, potentially allowing the space 
density to remain unchanged. From the measurements we 
have no constraints on how many of these galaxies there are 
and how they are distributed within the dark matter haloes 
and thus how they might change the clustering. Because the 
space density has changed in the merger model one could 
argue that we should compare our evolved high redshift 2pt 
correlation function with one measured from a sample of 
low redshift LRGs with a matching lower space density. The 
difficulty with this approach is deciding which galaxies to re- 
move from our observed sample in order to reduce the space 
density. 

An obvious choice would be to change the magnitude 
limit, thus removing the galaxies with the l owest stellar 
masses, equivalent to the approach taken in IWhite et al.l 



20071 ). However, in our merging model, we merge high red- 
shift central galaxies, and it seems unlikely that these would 
represent the LRGs with the lowest stellar masses. Alterna- 
tively, if we randomly sample the low redshift HOD, we will 
reduce the space density with out changing the clustering. 
This is equivalent to saying that the LRGs, which are newly 
formed by the merging of lower luminosity red galaxies at 
low redshift, trace the dark matter in the same way as the 
whole LRG population. If this is a true reflection of the evo- 
lution of the LRG population, then the randomly sampled 
measured HOD should look like the HOD produced by our 
central merging model. 

We show in Figure [T2] a comparison of the HOD of the 
best fit sub-Poisson merger model with the best HOD fit 
to the z = 0.19 measurement, along with the measured z 
= 0.19 HOD randomly sampled to match the space density 
of the merger model HOD. The left side of Figure [121 shows 
the HODs and the HODs weighted by the number density of 
the haloes in the same way as we've shown before. The right 
side shows the ratio of the HODs (top) and the difference 
between the weighted HODs (bottom) . For all but the lowest 
masses there is reasonable agreement between the randomly 
sampled HOD and the merger HOD. At the low mass end, 
the large difference is due in part to our having to force the 
z = 0.19 HOD to have a particular functional form; a form 
which the central merger model is not required to satisfy. 
There is still some discrepancy beyond that caused by the 
steps introduced by the binomial form of the sub-Poisson 
model, suggesting that any newly formed LRGs, which have 
been added to the low redshift sample, do not trace the dark 
matter in exactly the same way as the existing LRGs. 



7 COMPARISON WITH PREVIOUS WORK 
7.1 Merger Rates 

A number of authors have recently tried to co nstrain the 
merger rate of LRGs using a variety of methods. iBell et al] 
(2006) estimate that 50% of massive galaxies (> 5 x 
IO^^Mq) have experienced a major merger since z = 0.8. 
They also show that the merger rate increases with red- 
shift and provide a fitting formula for this increase. Apply- 
ing this formula to the redshift interval we are consider- 
ing here yields a merger rate of 21 % between z=0. 55 and 
z=0.19. The merger rate defined bv lBell et al] l|2006l ) is the 
equivalent of th e change in spac e density we measure i.e. 
9.2%. However, iBell et all (|2006! ) sample has a space den- 
sity of 33x 10~*/i~^Mpc^ which is 20 times higher than ours 
and thus consists of galaxies with typically much less stellar 
mass. The merger rate is believed to increase with decreas- 
ing stellar mass so any direct comparisons between the two 
mea surements are difficu lt . 

iMasiedi et al.l l|2006l ) use the small scale clustering to es- 
timate an LRG-LRG merger rate of 0.625% Gyr"^ for SDSS 
LRGs at, z = 0.25. This would correspond to 2% from z — 
0.55 to z — 0.19 far lower than our measurement. Applying 
the fitting f ormu la for the evolution of the me rger rate from 
Bell et al.l (|2006l ) normalised to match the iMasiedi et al.l 
(|200l ) value at 2 = 0.25 yields a rate ~4%, still a factor of 
2.5 lower than our best fitting value. Once again the galaxy 
samples aren't directly comparable since the space density of 
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Figure 12. In each panel the blue solid line shows the best fit to the 2SLAQ selection matched sample at z = 0.19. The effect of applying 
the sub-Poisson central merger model to the z = 0.55 HOD is shown as the red solid line and the z = 0.19 HOD fit random sampled 
to match the space density of the merger model is shown as the dashed blue line. The mean number of LRGs per halo as a function of 
halo mass is shown top left and the mean number of LRGs per halo times the number density of haloes as a function of mass is shown 
bottom left. The ratio of the numbers of LRGs is shown upper right and the difference in the space densities shown bottom right. 



LRGs in the iMasiedi et all l|200d ) sample are a factor of 3.5 
lower than the sample we use here, so one wo uld expect the 
merger rate to be lower for the more massive iMasiedi et al.l 
l|2006l) LRGs. 

IConrov. Ho. fc Whit3 (|2007l ) use N-body simulations to 
follow the accretion of halos sufficiently massive to host 
LRGs. They then compare this accretion history with the 
observed multiplicity function of LRGs at « ~ 0.3 (|Ho et al.l 
in order to constrain the LRG merger timescale and 
hence merger rate. Using this method they find a LRG-LRG 
merger rate approximately a fact or of two higher than that 
measured bv lMasiedi et al.1 l|2006l ') using the small scale clus- 
tering of LRGs. 

Using a ver y similar methodology to our own, 
I White et all (120071 ) estimate that ~ 1/3 of the low-redshift 
satellite galaxies must be destroyed (e.g. merge) in order to 
match the clustering evolution of luminous red galaxies be- 
tween z = 0.9 and z = 0.5 in the NDWFS. This corresponds 
to a merger rate of 3.4% Gyr^^, which would be 10.6% over 
our redshift interval. Th is rate is comprab le to our estimates; 
however, based on the iBell et al.l (|2006l ) trend, one would 
expect a factor of 2 i ncrease in the mea n rate due to the 
higher redshift of the IWhite et al.l I^OOT*) sample and also 
an increase due to the factor of 6 higher space density of 
their LRG s. There is, howeve r, one important difference be- 
tween the IWhite et al.l (|2007t ) study and the one presented 
here that may rectify some of the descrep e ncy in the merger 
rates. As mentioned above IWhite et al.l (|2003) adjust the 



space density of the low redshift LRG sample HOD fit with 
which they compare to their evolved high redshift sample. 
This is accomplished by adjusting the mass scale of the HOD 
fit by 7% to higher masses. This approach, of course, would 
reduce the space density and increase the clustering, result- 
ing in a lower amount of merging required to reduce the 
clustering produced by the passive evolution model to the 
measured level. Reducing the fraction of high-redshift cen- 
trals allowed to merge in the model similarly increases the 
clustering but also decreases the space density. 

Therefore there is only one unique combination of 
mass scale shift and merger rate that will match both the 
clustering and space density simultaneously. We find that 
increasing the z = 0.19 HOD mass scale by 6% and allowing 
63% of the high-redshift centrals to merge yields a large- 
scale bias of 1.93 and space density of 1.52x 10~*/i~^Mpc^ 
for both the measured low-redshift HOD and the evolved 
high-redshift HOD. This corresponds to a merger rate 
of 7.5±2.3% between z = 0.55 and z = 0.19. Figures [13] 
and [2] show the HOD and clustering respectively. Within 
the errors the merger model yields a good match with 
the measured HOD although the small scale clustering is 
a slightly poorer fit than the model with more merging 
shown in Figure 1111 This value is now in better agreement 
to that which one mi ght derive from the measurement 
of Masiedi et al.l (|2006l ) and the estimate of IWhite et al.l 
(|2007l ) although it still seems marginally higher. This may 
of course be due to the uncertainty in the dependence of 
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Figure 13. The mean number of LRGs per halo as a function of 
halo mass at z = 0.19 (black) for the 2SLAQ selection matched 
sample with the mass scale increased by 6.4%. The effect of pas- 
sively evolving the z = 0.55 fit to z = 0.19 is shown as the red 
line and the effect of including merging of the central galaxies is 
shown as the green line. 
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Figure 14. The ratio of the projected correlation functions to 
the best fit HOD at z = 0.19 where the mass scale of the HOD 
has been increased by 6.4%. The effect of passively evolving the 
z = 0.55 fit to z = 0.19 is shown as the red line and the effect 
of including merging of the central galaxies is shown as the green 
line. 

the merger rate with redshift and mass. Al t ernati vely the 
possible discrepancy with the IWhite et all (|2007[ l result, 
which uses a very similar method coupled to N-Body 
simulations could point to a deficiency in the current the- 
oretical models of the conditional mass function used herein. 

Finally, iMcIntosh et al.1 (|2007l 'l search for evidence of 
disturbance in close pairs of massive galaxies in z < 0.12 



groups to estimate the merger rate. They find that most 
of the mergers are occurring between approximately equal 
mass red progenitors and typically involve the central group 
galaxy, a picture that is consistent with the model we present 
here. T hey determine a mer ger rate 2-9 times higher than 
that of iMasiedi et aD (|2006l l for comparable galaxies and 
suggest that this is because their minimum group mass is 
3.5 X 10^'' Mq, higher than the typical halo mass of LRGs, 
and therefore the merger rate of LRGs increases with in- 
creasing halo mass. We show in Figure [15] the merger rate 
as a function of halo mass for our three merger models. This 
figure does indeed indicate a rapid increase in the merger 
rate in haloes with mass up to 3 or 4x10^''A/q, but with a 
decrease at higher masses. 

7.2 Semi-analytic Models 

lAlmeida et al.l (|2008l ) present a comparison of semi-analytic 
galaxy formation models to various properties of samples 
of LRGs very similar to the ones presented he re. They find 
that the one of their models ([Bowe r ot al. 2006) gives a good 
match to the luminosity function of SDSS LRGs ai z — 0.24, 
but over predicts the a bund ance of 2SLAQ LRGs at z = 
0.55. The iBower et all (|2006l ) model is also able to repro- 
duce the clustering of samples at both z = 0.5 and at z = 
0.24. They also present HODs generated from their models 
and compare them to the best fitting HODs for the samples 
presented within this paper, where the HOD fits are made 
using the same cosmological parameters as are used in the 
semi-a nalytic models. We only consider here the Bo wer et al.l 
(|2006l ) z=0.24 HOD, which is shown in Figure [16] since it 
matches both the LF and clustering of SDSS selected LRGs 
at z — 0.24. The plotted HOD has a quite difTerent form 
from the one we measure and is not reproducible with the 
formulation we have used in this paper. 

In addition these models predict satellite fractions of 
20-30% which is a factor of 2-3 times higher than our HOD 
fits yield, but a merger rate for the 2SLAQ selected sample 
of ~5% over our redshift range, in good agreement with the 
observations. Plotting the centra l and satellite HODs sepa- 
rately for the lBower et al.l l|2006l ) model (Figure [16)) demon- 
strates the reason for the high satellite fraction. There are 
many haloes that do not have an LRG central but do have 
satellite LRGs. It may appear surprising that the central 
galaxy within a halo does not meet the LRG selection cri- 
teria. Although the central galaxy is the most massive in 
terms of stellar mass and cold gas mass, it is not necessarily 
the brightest in the observer frame r-band. 

A more likely scenario, however, is the case in which 
the central galaxy is the brightest galaxy in t he halo, but 
does n ot match the LRG colour selection. In the lBower et al.l 
(|2006l ) model the suppression of gas cooling by AGN heating 
ramps up gradually from intermediate mass haloes, so some 
gas is still cooling in haloes with Mhaio ~ 10^^ Mq and 
being directed onto the central galaxy. This supply of cold 
gas results in recent star formation in the central galaxy. 
In more massive haloes, the cooling fiow is suppressed more 
strongly, so central galaxies in these haloes experience no 
recent star formation. It may also be the case that there 
are too many red satellites due to the instant aneous strip- 
ping of the gas a galaxy experiences in the iBower et al.l 
(j2006l ) model. However, one does need to remain cautious 
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Figure 15. The LRG merger rate as a function of halo mass. 

with these compar isons since even though at z = 0.24 the 
iBower et ahl l|2006l ) model does match both the LF and clus- 
tering, it is unable to reproduce the evolution of the LF, 
suggesting that it is still lacking in some areas. Even so, it 
does suggest that the form of the HOD we are using may 
be too simplistic when a colour selection is included along 
with a luminosity cut. We will investigate this further in a 
forthcoming paper, w hich includes bot h a better treatment 
of the gas stripping (|Font et al.ll2008l ) and a refined AGN 
feedback model. 
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Figure 16. The mean number of LRGs per halo as a function of 
halo mass (top) and the mean number of LRGs per halo times the 
number density of haloes as a function of mass (bottom) f or the 
Bower et al. model of LRGs presented in lAlmeida et al ] l|2008l 'l 
(red lines) compared to those generated from the fits to correla- 
tion functions present herein (blue lines). The total, central and 
satellite contributions are shown by the solid, dashed and dotted 
lines respectively. 



8 SUMMARY AND CONCLUSIONS 

We present here a detailed analysis of the clustering of Lu - 
min ous Red Galax ie s (as defined bv lEisenstein et al.l ()200ll ) 
and ICannon et al.l (|2006D ') as a function of redshift using 
samples of LRGs matched to have the same intrinsic colours 
and luminosities assuming passive evolution of their stellar 
populations. These galaxies represent the most massive in 
the universe with stellar masses lager than 10^^ Mq and 
space densities of ~ 10~"'/i^Mpc~^. We find that: 

• The amplitude of the clustering (ro) does not signif- 
icantly evolve with redshift over 0.15 < z < 0.6, whereas 
there is a marginally significant decrease in the slope (7) 
with increasing redshift. 

• The lack of evolution in the clustering amplitude on 
large-scales is inconsistent with a picture in which the 
LRGs have purely passive evolution undergoing no major 
mergers over this time period, and rules out this passive 
model at 98% significance. 

• A IfOD where the fraction of haloes which host 
central galaxies {Nc\A'I) = exp(— Afmin/M) and only 
haloes which host centrals can host satellites where the 
satellites are drawn from a Poisson distribution with mean 
{Ns\M) = (M/Mi)" is able to accurately reproduce the 
clustering and space density of our LRG samples. Within 
this framework the LRGs are predicted to be hosted in 
haloes with a typical mass close to 10^"^ Mq which 



increases by ~ 50% from z — 0.55 to z — 0.2, and to have 
satellite fractions increasing from ~ 5 to 10% over this 
time. The LRGs are found to have a bias ~ 2 and which 
decreases with redshift at a much greater rate than would 
be predicted for the passive no merger case. 

• We introduce an analytic approach to describe the 
evolution of the HOD with redshift, and demonstrate 
that this guarantees that the comoving density remains 
constant and the large scale bias evolves in accordance with 
the continuity equation. We use this approach to further 
demonstrate that the passive evolution of the LRG HOD 
from z = 0.55 is inconsistent with the measurements at z 
= 0.19 at greater than 99.9% significance, predicting far 
too many satellite galaxies at z — 0.19 and greatly over 
estimating the clustering strength on all scales. 

• We introduce a model in which high-redshift centrals 
are allowed to merge with other high-redshift centrals 
occupying the same halo at low-redshift. This choice is 
motivated by the fact that centrals are likely to be more 
massive than satellites, so dynamical friction may be 
more effective at making these objects merge with the 
true low-redshift central. This model is able to accurately 
match the large-scale clustering evolution of the LRGs. We 
demonstrate that the small scale clustering is dependent on 
the parametrisation of the scatter in halo merger histories. 
We investigate three models for this scatter and find that 
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both the sub-Poisson and Poisson models are able to match 
the small-scale clustering evolution. However, the Poisson 
model requires a much larger LRG-LRG merger rate (20%) 
which is not favour ed by either the evolution of the LRG 
luminoisty function l|Wake et al.ll2006l ) or other independent 
measures of the LR G-LRG merger rate l|Masiedi et al.|[2006l : 
IWhite et al.]|2007l ). We therefore favour the best motivated 
sub-Poisson scatter giving observational support to this 
generic prediction of hierarchical models. 

• In order to match the clustering evolution we require 
an LRG-LRG merger rate of 7.5±2.3% from z = 0.55 to z = 
0.19 corresponding to 2.4% Gyr~^. This is probably consis- 
tent with other measurements of the merger rate of massive 
red galaxies given the uncertainties in how the merger rate 
depends on the mass of the galaxy and evolves with redshift. 

• Although some merging is required to match the 
clustering evolution, the merger rate is sufficiently small 
that it is entirely compatible with the low rate of evolution 
in the luminosity function of LRGs found in IWake et al.l 
l|2006t ). 

• We compare in detail the measured HOD for one of the 
LRG samples to that predicted by the latest semi-analytic 
models of galaxy formation for a very similar sample 

of LRGs which matches bot h the luminos ity function 

and clustering as described in lAlmeida et al.1 (|2008l ). The 
model HOD is very different to our fit, and would not 
be reproducible by the functional form of the HOD we 
assume. In particular the model has many haloes that 
contain LRG satellites where the central is not an LRG. 
This suggests that a more sophisticated form of the HOD 
may be required for galaxy samples selected by colour in 
addition to luminosity, although caution is required as the 
semi-analytic model is still unable to accurately reproduce 
the evolution of the LRG population. 

• Our halo model analysis of the relation between the 
low- and hi gh-redshift p opulations is similar i n spiri t 
to those of IWhite et all (|20o3) and ISeo et al.l l|200it ). 
However, whereas their work used numerical simulations, 
our approach is entirely analytic. This means that our 
analysis relies heavily on the accuracy of current models 
of N{m\M), the conditional mass function. These models 
are not particularly ac curate for small redshift intervals 
ISheth fc TormenI ||2002D . so we hope that our analysis will 
generate interest in improving these models. 

• Our analysis also highlights the need for a better 
understanding of the stochasticity in halo merger histories. 
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APPENDIX A: CONSTANT COMOVING 
NUMBER DENSITY IN THE HALO MODEL 

Al Large scale clustering in real-space 

Let g{m) denote the mean number of galaxies in haloes of 
mass m at some early time, and let G{M) denote a similar 
quantity at some later time. If haloes merge but galaxies do 
not, then 



G{M)= / dmN{m\M)g{m) 



(Al) 



where N{m\M) denotes the mean number of m haloes from 
the earlier epoch which have been incorporated into M 
haloes by the later epoch. Lacey & Cole (1993) and Sheth & 
Tormen (2002) discuss models for N{m\M) that are consis- 
tent with the halo abundances of Press & Schechter (1974) 
and Sheth & Tormen (1999). 

To see that the number density of galaxies has indeed 
not changed, note that 



dMn{M) G{M) 



AMn{M) I dmN{m\M)g{m) 
Jo 

dmg(m) / dM n(M) N{m\M) 



dmn(m) g(m). 



(A2) 



The first equality expresses the number density as an inte- 
gral over the low-redshift halo population, whereas the final 
equality integrates over the high-redshift population. The 
associated large scale bias factor at the later time is 



6o - 1 = 



"dM<MMMly,^M)-l] 



dm 



g{m)n(m) 



, , ^ , n(M) N (m\M) „ , ^ ,^ , , 
X / dM-i — \ V ' — '-\b(M)-l]. A3) 
n(m) 



Now, 

^ , _ dlnn(M) 
dOc 



(A4) 



(Sheth & Tormen 1999) and the algebra in Abbas & Sheth 
(2005) shows that the expression above reduces to 



ho - 1 



dm 



g{m)n(m) [fe(m) — 1] 



[h, - 1)/{Dq/D,) 



(A5) 
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where D is the linear theory growth factor. (If the later 
time is the present in an Einstein de-Sitter universe, then 
Dq/Dz = ao/oz = 1 + 2.) This shows explicitly that the 
halo model calculation of the evolution of the bias in the no- 
mcrgcr model is the same as that derived from an argument 
based on the continuity equation (Nusser & Davis 1994; Fry 
1996). Note that the bias factor evolves even though the 
number density does not. 

One might wonder if, although the bias factor evolves, 
the clustering strength itself docs not. The ratio of the large- 
scale clustering signal at the two epochs is 

e„(r)_bgg§_/ 6o V 

Ur) blDl \bo-l + D,/DoJ ' ^ ' 

since < Dg, the later epoch is more strongly clustered. 

For example, for 60 — 2 and Dz/Do = 2/3, this factor is 
(6/5)^ = 1.44. Setting D^/Do < 1 illustrates a fact that is 
often overlooked: the clustering strength of highly biased ob- 
jects (i.e., the most massive haloes) evolves very little, even 
though the clustering of the dark matter itself has evolved 
significantly: (Dq/D^)^ ^ 1. The most massive objects do 
not move far from their initial comoving positions. 

This calculation suggests a simple test of the null hy- 
pothesis that two populations having the same comoving 
number density are related by the no-merger evolution 
model: If the measured clustering signal has not evolved, 
or if the high redshift sample is more strongly clustered, 
then the hypothesis can be rejected. 



A2 Small scale clustering in real-space 

The continuity equation argument is restricted to the large 
scales on which linear theory applies. The virtue of writing 
this in terms of halo abundances is that it shows clearly how 
to extend the model to predict the clustering signal in the 
no-merger model even on small scales. In particular, two ad- 
ditional pieces of information are required: a model for how 
the galax;ies are distributed around the centre of their parent 
haloes, and the second factorial moment G2{M) of the dis- 
tribution p{N\M) of the number of galaxies A^^ at fixed halo 
mass M. Sheth et al. (2001) show that, on scales larger than 
approximately half a Megaparsec, it is more important to 
model the first two moments Gi(M) and G2(M) accurately 
than the density profiles; in particular, the approximation 
that the spatial distribution of the galaxies is the same as 
that of the dark matter is sufficiently accurate. Hence, if 
we know the second factorial moment of how galaxies popu- 
late haloes, then we can describe the no-merger correlation 
function on small scales as well. Simulations indicate that 
in haloes which host more than one galaxy, p{N — 1|M) 
is a Poisson distribution with mean Gi — 1. This specifies 
G2(M). 



